ELSEVIER 


Science of the Total Environment 


Science of the Total Environment xxx (2017) xxx-xxx 


Contents lists available at ScienceDirect 


journal homepage: www.elsevier.com 


Seagrass blue carbon dynamics in the Gulf of Mexico: Stocks, losses from 
anthropogenic disturbance, and gains through seagrass restoration 


Anitra Thorhaug” *, Helen M. Poulos” °, Jorge Lopez-Portillo’, Timothy C.W. Ku‘, Graeme P. Berlyn® 


* Yale School of Forestry and Environmental Studies, New Haven, CT 06511, United States 

> Wesleyan University, College of the Environment, Middletown, CT 06457, United States 

© Instituto de Ecologia, A.C. (INECOL), Xalapa, Veracruz 91070, Mexico 

: Wesleyan University, Department of Earth and Environmental Science, Middletown, CT 06459, United States 


ARTICLE INFO 


Article history: 

Received 17 April 2017 

Received in revised form 21 June 2017 
Accepted 22 June 2017 

Available online xxx 


Editor: Jay Gan 


Keywords: 

Blue carbon 

Seagrass carbon 

Seagrass restoration 

Gulf of Mexico 

Seagrass impact carbon loss 


ABSTRACT 


Seagrasses comprise a substantive North American and Caribbean Sea blue carbon sink. Yet fine-scale es- 
timates of seagrass carbon stocks, fluxes from anthropogenic disturbances, and potential gains in sedimen- 
tary carbon from seagrass restoration are lacking for most of the Western Hemisphere. To begin to fill this 
knowledge gap in the subtropics and tropics, we quantified organic carbon (C,,.) stocks, losses, and gains 
from restorations at 8 previously-disturbed seagrass sites around the Gulf of Mexico (GoM) (n= 128 cores). 
Mean natural seagrass C,,, stocks were 25.7 + 6.7 Mg C,,, ha | ' around the GoM, while mean Corg Stocks at 
adjacent barren sites that had previously hosted seagrass were 17.8 Mg Co,. ha” ' Restored seagrass beds con- 
tained a mean of 38.7 + 13.1 Mg Cy, ha” "Mean Core losses differed by anthropogenic impact type, but aver- 
aged 20.98 + 7.14 Mg C,,, ha” a Cor gains from seagrass restoration averaged 20.96 + 8.59 Mg ha ' These 
results, when combined with the similarity between natural and restored C,,, content, highlight the potential 
of seagrass restoration for mitigating seagrass C,,, losses from prior impact events. Our GoM basin-wide es- 
timates of natural C,,, totaled ~ 36.4 Tg for the 947,327 ha for the USA-GoM. Including Mexico, the total 
basin contained an estimated 37.2-37.5 Tg Cy.9. Regional US-GoM losses totaled 21.69 Tg Cory. Cory losses 
differed significantly among anthropogenic impacts. Yet, seagrass restoration appears to be an important cli- 
mate change mitigation strategy that could be implemented elsewhere throughout the tropics and subtropics. 


1. Introduction 


Seagrasses are a critical, but continually-diminishing blue carbon 
foundation habitat. The IPCC (Intergovernmental Panel on Climate 
Change) Wetlands Committee only recently recognized that wetlands 
and submerged aquatic vegetation including marshes, mangroves, and 
seagrasses contribute significantly to the global stored carbon sink, 
and Hiraishi et al. (2014) called upon governments to implement 
more sophisticated inventories of wetlands carbon emission factors 
to move beyond the use of basic default values for tidal wetlands 
(140 gy _') and basic equations (i.e., Tier I standards). Current base- 
line estimates indicate that global seagrass carbon pool lies between 
4.2 and 8.4 Pg of organic carbon (C,,) (Fourqurean et al., 2012a). 
Yet, these important seagrass blue carbon sinks continue to be lost at 
an alarming rate (Waycott et al., 2009). Given that blue carbon se- 
questration rates can be orders of magnitude higher than terrestrial 
carbon pools (Fourqurean et al., 2012a), information about seagrass 
blue carbon cycling is an important factor in IPCC decision-mak- 
ing and carbon accounting. Presently, regional subtropical and 
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tropical seagrass blue carbon stock data remain sparse for seagrass and 
we lack a coherent fingerprint of the spatial variation in C,,, stocks 
in relation to local environmental influences. Nor do we understand 
the potential of seagrass restoration for bolstering blue carbon stocks 
in these regions. Such information is critical for achieving long-term 
Tier II (regional/country-specific emission factor values) and Tier III 
(more complex, site-specific emission factor values) blue carbon ac- 
counting goals and REDD + (reducing emissions from deforestation 
and forest degradation) carbon offset eligibility. 


1.1. Seagrass blue carbon losses from anthropogenic disturbances 


Cor, fluxes from seagrass die-off after anthropogenic disturbance 
events can vary considerably depending on local environmental con- 
ditions, C,,,. inputs to the estuary, and the type and magnitude of the 
disturbance (Macreadie et al., 2013). Pendleton et al. (2012) gener- 
ated a rough global estimate of blue carbon emissions from conver- 
sion and degradation of coastal marine vegetation. However, their es- 
timates utilized default mean C,,, stock values of 140 Mg ha | from 
Fourqurean et al. (2012a), and they did not incorporate spatial varia- 
tion in sedimentary carbon stocks or the effects of different types of 
impacts on seagrass C,,,-cycling. Several researchers have presented 
site-specific data on seagrass carbon flux from human impacts includ- 
ing Marba et al. (2015) in Western Australia on the Indian Ocean 
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at a location where the seagrass ecosystem sediment collapsed result- 
ing in the disappearance of the seagrass carbon and its associated sed- 
iments in response to disturbance. Likewise, Macreadie et al. (2014) 
experimentally disturbed the sediment in small plots, and measured 
sediment carbon content over 24 months in central eastern Australia, 
but found no significant Co,, loss. Such results highlight the need 
for more studies that address the mechanisms of seagrass carbon cy- 
cling following anthropogenic disturbances. The few existing empir- 
ical studies on seagrass blue carbon losses suggest that making com- 
mon assumptions about coastal marine carbon cycling can be mislead- 
ing (Macreadie et al., 2014). 


1.2. Seagrass restoration as a climate change mitigation strategy 


van Katwijk et al. (2016) demonstrated the importance of seagrass 
restoration on a global-scale for restoring coastal marine ecosystem 
services, yet the effects of restoration activities on blue carbon cycling 
have yet to be examined fully. Seagrass restoration theoretically pro- 
vides one potential solution for mitigating the losses of blue carbon 
stocks from anthropogenic disturbances. Yet, few prior studies have 
examined the potential of seagrass restoration for bolstering blue car- 
bon stocks. To date, data exist from only two sites (Chesapeake Bay, 
USA (Greiner et al., 2013) and Perth, Western Australia (Marba et al., 
2015)) that quantify the effects of seagrass restoration on seagrass sed- 
imentary carbon stocks. The results of both studies indicate that C,,. 
stocks in restored areas following seagrass die off from anthropogenic 
impacts can vary considerably. These studies represent two disparate 
data points from one temperate high energy site on Western Australia's 
Indian Ocean and one temperate Atlantic estuarine site in Chesapeake 
Bay, USA. Regional-scale studies of the effects of multiple seagrass 
restorations on blue carbon stocks under a variety of environmental 
and anthropogenic stressors are clearly needed for achieving detailed 
and regional-scale Tier II and III IPCC carbon accounting goals. 

This study builds upon the limited body of literature on seagrass 
sedimentary carbon accounting via a regional-scale analysis of sea- 
grass Co,, stocks, losses from an array of anthropogenic disturbances, 
and C,,, gains through multiple seagrass restorations across the USA 
GoM. Our overarching goal was to examine the dynamics of seagrass 
sedimentary carbon flux from disturbance and to evaluate the potential 
of subtropical/tropical seagrass restoration as a climate change mit- 
igation strategy. We hypothesized that the effects of different types 
of seagrass disturbances on C,,. losses would vary in accord with 
the disturbance type, and that different types of disturbances would 
also influence the subsequent carbon sequestration potential of sea- 
grass beds that were then restored after such impacts. We also hy- 
pothesized that seagrass restoration would bolster blue carbon stocks 
in barren areas where seagrass die-off had previously occurred. We 
tested these hypotheses by systematically collecting C,,, data from 
across 8 subtropical/tropical seagrass sites that had all been restored 
using a consistent restoration methodology across the USA GoM, 
which is considered a of North American carbon hotspot (Hofmann 
et al., 2008; Herrmann et al., 2015). Each sample site had experi- 
enced major seagrass die-off from an anthropogenic disturbance, and 
then Thorhaug subsequently restored portions of each site using stan- 
dardized restoration methods over the last 4 decades in a series of 
restorations (see Table | for descriptions of individual restoration pro- 
jects). The present study uses known seagrass kill dates (1.e. the dates 
of closing of thermal and sewage effluents, vessel groundings, and 
artificial-fill dates of US Army Corps of Engineers) and the exact 
restoration dates to estimate Co, losses and gains at 8 sites across the 
USA GoM (Table 1 contains descriptions of each study site and its 


Table 1 
Site information for eight Gulf of Mexico seagrass sample sites. All sites contained the 
four treatment blocks including natural seagrass, restored seagrass, impacted barren ar- 
eas that were never restored, and always barren sites. The age (y) of each restoration, 
species composition, spatial location, sediment type, and impact type are listed for each 
location. 


Impact type with 

Restoration Sediment age; restoration age 

age & species Sites & location type, GIS & citations 

42 y Thalassia Turkey Point, S. Biscayne Medium Heated effluent 
testudinum, Bay, FL. Abundant quartz sand = (1969-71) Bader et 
Halodule Mangroves al. (1970). Planting 
wrightii 252,510.88 N 801,934.09 W 1973 by Thorhaug 

(1974). 

32 y Thalassia 24—27th Street NEN. Course Sewage (ended 1958) 
testudinum, Biscayne Bay, FL quartz sand = =McNulty (1970) 
Halodule 254,352.35 N 801,233.19 W Planting 1983, 
wrightii, Thorhaug (1985) 
Syringodium 
filiforme 

37 y Thalassia N. Julia Tuttle (36th St) Rock and Dredge & filled 
testudinum, Causeway, N. Biscayne rubble construction (1959) 
Halodule Bay, FL Mangroves present _ overlain Chardon (1976). 
wrightii, 254,840.54 N 801,003.66 W with quartz Planting 1978 by 
Syringodium sand. Thorhaug (1980a) 
filiforme Constructed 

causeway. 
29 y Thalassia Grove Isle, Central Biscayne Quartz sand Dredge and fill (early 


testudinum, Bay, FL, Sparse mangroves 1961), Harlem 
Halodule present (1979), Planting 
wrightii, 254,342.11 N 801,302.72 W 1985 by Thorhaug 
Syringodium (1987) 

filiforme 


15 y Halodule Predator Island, E. Northern Quartz Sand. Fill for 


wrightii, Laguna Madre, TX Scrape Channelization, 

Halophila Adjacent marshes. down of (1947) artificial 

decipiens 273,220.22 N 971,113.41 W former island scrape-down 
dredge (1998). Planting 
island. 1999, Thorhaug 


(2001) 


1 y Halodule — Shamrock Shoreline, East Sand, quartz, Filled causeway & 


wrightii Corpus Christi Bay, TX. oyster shell — shoreline & oil spill 
Adjacent marshes. plus clay. 1942's. Planting 
274,516.74 N 970,927.16 W 2013 by Thorhaug 
(2014) 
2 y Halodule San Luis Pass, Galveston Sand. Fill 1969 (TX 
wrightii Bay, TX. Adjacent marshes. Quartz. Handbook) 
290,551.61 N 950,638.06 W Planting, 2013 
Thorhaug (2014) 
ly Halodule —_Aransas/St Charles Bay, TX. Sand quartz, Lg. Salinity 
wrightii 280,851.94 N 965,819.81 W with marsh (2011-3), bridge 
Peat construction 
underlying (2011-13), 
Shoreline fill 


(2005-2013). 
Aluminum railings. 
Planting 2013 by 
Thorhaug (2014) 


anthropogenic impact). The dates of both the impact event and the 
restoration date, provided a fine-scale clock for quantifying both Co, 
losses since disturbance and gains since restoration. 


2. Study region 
2.1. Gulf of Mexico seagrasses 


The GoM is an internal North American basin, that is dominated 
by sedimentary estuaries that contain seagrasses as a foundation habi- 
tat throughout the basin. Soft sandy and muddy bottoms are char- 
acteristic of estuaries across the region, and these wetlands support 
large expanses of seagrasses, marshes, and mangroves. Two ben- 
thic western Florida shelves also support seagrass. Seagrasses extents 
in the USA GoM have been estimated at near | million ha (chiefly 
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found in Florida and in Texas's Laguna Madre (Appendix A) (Duke 
and Kruczynski, 1992; Pulich and Hardegree, 2006; Handley et al., 
2007; Yarbro and Carlson, 2011) with rough estimates of 25,000 ha 
in the Mexican portion (Green and Short, 2003). Seagrass productiv- 
ity is high in the GoM due to its year-long growing season and clear 
waters (Thorhaug et al., 1978; Dunton, 1994; Livingston et al., 1998). 
Extensive anthropogenic development has occurred chiefly during the 
last century, which has decimated intertidal wetlands and seagrasses 
throughout the region. Duke and Kruczynski (1992) estimated that 
54% (19,250 km?) of all wetlands and seagrasses in the GoM were lost 
by 1992. To date, this seagrass loss was considered to chiefly affect 
fisheries and habitat extent. However, such losses likely also have ef- 
fects on GoM seagrass carbon capture and storage. 


3. Methods 
3.1. Field methods 


Field sampling took place at low tide during the summer and early 
fall of 2014 at all sites. At each sample site, quadruplicate 7-cm ra- 
dius and 20-cm deep sediment cores were collected from the following 
treatments: (1) “always barren” (AB) over 50 y; (2) naturally-occur- 
ring seagrass (NS) for over at least 50 y; (3) “impacted barren” (IB) 
former seagrass meadows killed by the disturbance but receiving no 
restoration and is still barren of seagrass in 2014, and (4) restored sea- 
grass (RS) (Table 1). The 16 cores at each sample site (4 per treat- 
ment) were taken at 7.5-m intervals along a transect that ran through 
the center of each treatment type and parallel to the current at ~ 0.75 
to 1m water depth (n= 128 cores). Transects were run through the 
center of each treatment type to minimize edge effects from prior sea- 
grass die-off. All four treatments (AB, NS, IB, and RS) were present 
in patches within each sample site, and the NS, IB, and AB patches at 
each site were all located within a 650-m radius from each restoration 
site. Temporal changes in seagrass cover were estimated prior to go- 
ing in the field using historical aerial imagery (https://coast.noaa.gov/ 
dataregistry/search/collection). Seagrass die-off, or loss, generally oc- 
curred simultaneously at each site at a known date in response to a 
known disturbance event. These seagrass die-off dates or “clocks” in- 
dicated the time of the vessel (2 sites), fill (3 sites), stoppage of efflu- 
ents (2 sites), or salinity modification event (1 site) that impacted each 
seagrass meadow, after which natural sediment accumulation presum- 
ably resumed. A second “clock” was the known time of planting each 
site. 

We chose a coring depth of 20 cm based on the rooting depths of 
Thalassia testudinum and Halodule wrightii and published accretion 
rates for Biscayne Bay, Corpus Christi, and Galveston Bay in nearby 
seagrass beds (Wingard et al., 2003; Yeager et al., 2006; Ravens et al., 
2009), which suggested that a 20 cm deep core would capture the car- 
bon lost from a 60-year-old impact event and the carbon gained from a 
40-year-old restoration (the oldest restoration). Our field observations 
further justified this sediment core depth because the anthropogenic 
impact and subsequent sediment restoration layers were visible in our 
restored seagrass sediment cores as large changes in sediment texture, 
composition, and color at various sample depths at sites within the 
20 cm IB and RS cores. 

Samples were immediately extruded after coring, put on ice in the 
field, and then immediately refrigerated at -4 °C upon return from the 
field until processing immediately after collection to prevent micro- 
bial processes from metabolizing the carbon. Sediments in most cores 
were aluminosilicate sand and they experienced little compaction from 
coring. However, some cores did show evidence of sediment com- 
paction (11% of all cores taken), and for those cores, we applied 


compaction correction factors prior to sectioning cores for further sed- 
iment analysis following Howard et al. (2014). Cores were sectioned 
and homogenized immediately upon return from the field at three sam- 
pling intervals per core including 0-5 cm, 10-15 cm, and 10-20 cm 
depths (n= 384). Shells, rocks, and large rhizomes were removed 
prior to subsequent C analyses. 


3.2. Lab methods 


In the lab, samples were transferred into pre-weighed 15 mL cen- 
trifuge tubes, and wet sediment volumes were recorded after centrifu- 
gation at 3000 rpm for 8 min. Samples were then dried at 70 °C to a 
constant mass and the dry bulk density (BD) was estimated by divid- 
ing the dry mass of each sediment section by the wet volume. Dried 
0-5 cm, 5—10 cm, and 10-20 cm samples were homogenized using a 
mortar and pestle and transferred to 20 mL glass scintillation vials 
(n = 384). Carbon analyses followed procedures outlined by Howard 
et al. (2014). Loss on ignition (% OM), and total and inorganic carbon 
(% Cio, % C;) were measured on a CE Elantech Flash 1112 Elemental 
Analyzer (n = 768). Inorganic carbon measurements were completed 
after a loss on ignition (LOI) measurement where ~ 500 mg subsam- 
ple was ashed in a muffle furnace at 500 °C for 5 h. Inorganic carbon 
results were then scaled to the original dry sediment masses by using 
ash to sediment weight ratios. The percentage of organic carbon (C,,.) 
was calculated as the difference between Ci and C,,;. Total C,,.. den- 
sity (Mg Cog ha” ) and total organic carbon standing stock were esti- 
mated using the methods described in Howard et al. (2014). 


3.3. Statistical analysis 


Data were examined for normality and transformed as necessary 
prior to analysis. We used mixed models and the nlme package of R 
for analyzing differences in % C and C,,, content among treatments, 
disturbance types, sample site regions, and restored seagrass bed ages. 
Differences in % Co, were also evaluated over the depth profiles of 
each treatment via mixed models. In all analyses, cores were nested 
within each site-treatment combination to account for the covariance 
structure in the nested sampling design. Significant differences were 
assessed among the following site and treatment combinations: 


Treatment-level differences in Corg (NS, RS, IB, vs. AB). 

Cor. depth by treatment differences in carbon content. 

Species-level differences in sedimentary C,,,, with core nested 
within treatment. 

Differences among sample sites in natural vs. restored C,,, by region 
and restoration age: young (< 3 years), intermediate (= 4-15 years), 
and old (> 15) since restoration. 

Differences in restoration response by impact or disturbance type 
(dredge and fill, effluent, turbidity, thermal pollution, or salinity). 


org? 


Sample sites and treatments were included as fixed effects. Ran- 
dom effects included the intercept, treatment-site relationships, and 
the core nested within sample site. The residuals of each model were 
inspected for deviations from homoscedasticity, and only models con- 
taining residuals without obvious deviations from normality were 
kept. We used linear regression to estimate the relationship between % 
OM and % C,,. to assess the effectiveness of the LOI technique. Dif- 
ferences in C,,, stocks between seagrass species (7. testudinum and H. 
wrightii) were tested via a Student's f-test. 

Losses in Co,g from anthropogenic impact and gains from restora- 


tion were calculated on a site-by-site basis using mean C,,, stock val- 
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ues for the four cores in each treatment at each sample site for the 
entire 0-20 cm core. We followed methods employed by Macreadie 
et al. (2013) to assess variation in sediment bulk density across the 
depth profiles of the cores in each treatment, which indicated no sig- 
nificant differences among sample depths or treatments (Appendix A), 
and supported the assumption that the carbon stocks to this depth were 
comparable among sites. Thus, we estimated the loss in Co, from prior 
impact events that resulted in seagrass die-off by comparing the four 
different treatments (AB, IB, NS, and RS) at each sample site. Loss in 
Corg stock in the top 20 cm of the sediment due to seagrass die-off at 
each site was estimated as: 


Loss = mean natural seagrass C,,, 


— (mean impacted barren Cy, 


— mean always barren Cae) ; 


This metric denotes the amount of carbon which has been lost via 
leakage or perturbation of sediment over the period between the im- 
pact when the seagrass bed was killed or removed by impact and our 
field sampling. The original natural seagrass bed at each site prior to 
the disturbance was the control for original bed and the spatial prox- 
imity of all four treatments within 650 m among each other makes this 
comparison possible. Our assumption here was that the impacted bar- 
ren sediments did not biologically-actively continue to sequester C af- 
ter seagrass was killed. 

We estimated the Co. gained through restoration as the mean Cor 
among cores at each site in restored beds minus mean C,,, in for 
the four impacted barren cores at each sample site. Values were re- 
ported as overall averages by treatment, impact type, and by restora- 
tion age (young - <3 years old, intermediate - 4-15 years old, and 
old - > 15 years since restoration). Estimates were calculated using 
site-specific, fine-scale estimates of total gain and gain per year using 
total Cy,. stocks for each 20-cm core and by taking the difference in 
Corg in restored sites from Co,, in natural seagrass sites at those same 
sites divided by the number of years since restoration. 


4. Results 


4.1. Seagrass organic carbon stocks among sites, restored vs. natural 
beds vs. barren 


Organic carbon stocks (C,,,) varied significantly among treatments 
and across the 8 USA GoM sites (F = 3.84, P = 0.011), and there was 
a highly significant relationship between LOI and % C,,. (Figs. 1-3; 
Table 2). The mean C,,, in the top 20-cm of the sediment for nat- 
ural seagrass beds at the sample sites was 25.7 + 6.7 Mg Cory ha - y 
Restored seagrass sites had higher mean Cy, — stocks 
(38.7 + 13.1 Mg ha 1) than natural beds, although differences were 
not significant. However, restored seagrass beds had significantly 
higher C,,, than impacted-barren and always-barren areas (P < 0.05). 
The baseline (control) of each estuary carbon stock, the “always bar- 
ren” samples, varied substantially among sites (from 4.11 to 
35.61 Mg Co ha ') and seagrass species (Fig. 2, Table 2), but the 
Corg stock sequestered in natural seagrasses was always above that 
barren estuarine baseline (AB) (Table 2). 


4.2. Changes in C 


org With sediment depth 


The organic carbon density by depth results (Fig. 4) indicated that 
substantial carbon was no longer present for the 0-15 cm of upper sed- 
iment layer depths in impacted barren (IB) areas compared to their 
former carbon levels in the natural seagrass (NS) beds. In the upper 
0-5 cm, IB Co, was not statistically different from AB. However, at 
5-10 and 10-20 cm depths IB C,,, was significantly higher than AB 
indicating that C,,, was no longer present for the 0-15 cm of the upper 
sediment layer depths in impacted areas IB but that residual Co,. was 
present in those cores at deeper sediment depths. 


4.3. Cog gains from restoration and restoration age effects 


Mean organic carbon gains through seagrass restoration were 
20.96 + 8.59 Mg Cog ha | across the sites for the top 20 cm of sedi- 
ment. We note variable values by sample site location and by the age 
of the restored bed (Fig. 5, Table 2). Co. in restored sites differed 
significantly by restoration age (young - <3 years old, intermedi- 
ate - 4-15 years old, and old - > 15 years since restoration) (F = 3.74, 


Total C.,, stocks (20 cm core) 


25 b b 


15 
10 
5 
0 
ab ib 


Mg C4, hat 
nN 
o 


a 


| 
ns rs 


Fig. 1. GoM seagrass C,,, content of natural seagrass (ns), restored seagrass (rs), always-barren (ab), and impacted-barren (ib) sediment. Values reported as means (+ SE) across the 
8 Gulf of Mexico estuarine sites. C,,, differs significantly among treatments per a mixed model analysis (P = 0.02). Pairwise comparisons are shown with letters. Treatments not 


sharing letters significantly different P < 0.05-levels. 
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® Naturalseagrass 


“ Impactedbarren * 


Fig. 2. Mean (+ SE) organic carbon content of 20 cm cores from H. wrightii and T. tes- 
tudinum natural seagrasses from 8 sampled Gulf of Mexico Sites vs. proximate previ- 
ously-impacted seagrass meadow now barren. Significant differences (P < 0.05) among 
treatments and species are depicted with asterisks (*) according to a t-test. 
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Fig. 3. Regression results between % C,,. of the loss corrected percent weight of sam- 
ples versus loss on ignition (LOI) from seagrass sediment samples from 8 sites in the 
Gulf of Mexico, including naturally occurring, restored seagrass, always barren and im- 
pacted barren sediment carbon values. 


P=0.033) and by location (F = 6.34, P= 0.016), but not by species 
(Figs. 2 and 5, Table 2). The intermediate age sites had the highest car- 
bon accumulation rates per y (Fig. 5, Table 2). 


4.4.C, 


org Losses from impact and disturbance 


Average organic carbon loss was 20.98 + 7.14 Mg Cog ha ! for 
the 8 sites across the USA Gulf of Mexico. Around the GoM basin, 
seagrass C,,. loss varied significantly among impact types, species, 
and location (Table 2, Fig. 6, Appendixes B and C). The highest 
Core losses (57.1 Mg Co. ha | ') occurred at the Corpus Christi Oil & 
Gas site vessel grounding in Northern Laguna Madre, Texas. The 
lowest losses were at the San Luis Pass, Galveston Bay, Texas 
(2.91 Mg C,,, ha '). 


org 


4.5. Pollutant-disturbance type effects 


Seagrass organic carbon loss also differed among impact/distur- 
bance types (Fig. 6). We also arranged differences among carbon into 
two groups: 1) where die-off from impact occurred once, but then 
seagrasses were conserved, and 2) sites that have continually experi- 
enced impact effects since initial seagrass die-off. The impact event 


type significantly influenced restored seagrass bed C,,, stocks 
(F = 21.2, P<0.0001). The site that experienced a sewage spill (36th 
St.) had the highest C,,, gain following seagrass loss, whereas the high 
salinity site had lowest C,,. stock following restoration (Aransas Bay). 
Seagrass meadows exposed to thermal pollution or vessel ground- 
ing, had intermediate C,,, stocks that were similar among sites, al- 
though the vessel grounding occurred in south Texas (colonized by H. 
wrightii) and thermal pollution occurred in South Florida (colonized 
by T. testudinum). Both sites are now preserves, surrounded by nat- 
ural shorelines and remain relatively undisturbed. The natural seagrass 
Corg at sites that were continually impacted even after the restorations 
had substantially lower C,,, stocks (mean 11.7 Mg C,,. ha ') than the 
once-impacted then-undisturbed sites (mean 42.5 Mg C,,, ha ') both 
in Texas and Florida (Turkey Point, Predator, and Grove Isle). Nat- 
ural seagrass C,,. stocks were substantially lower than the Co, stock 
of barren sediment at the only once-impacted sites (Table 2). 


4.6. Regional total C4. sequestration 


org * 


The mean natural seagrass C,,, sequestration was 
39.2 2.2 Mg Cy. ha” ' for Florida (total Cog = 34.7 Tg) and 
13.9 + 8.3 Mg Cory ha ' in Texas (total Corg 1.5 Tg). The average 
Texas Co, stock across sites that were restored 2 to 15 y ago was 
lower than older South Florida restored sites (35-42 y ago). Back- 
ground C,,, stocks in always barren sediments also differed signifi- 
cantly among disturbance-impact types and sites (Fig. 6, P < 0.05). 

Florida had the highest regional Co,, stock (21.6 Tg Cog) based on 
a seagrass extent of 850,874 ha (Appendix B), followed by Texas with 
94,389 ha of seagrass (2.5 Tg C,,.), and Alabama, Mississippi and 
Louisiana's combine extent of 1964 ha (0.049 Tg C,,.). The seagrass 
Corg stock for the entire USA GoM was estimated by multiplying the 
mean C,,. stock per ha by the federal/state government-measured spa- 
tial extents (see Table Appendix B for sources of spatial extent data) 
to obtain estimates of total regional C,,, stock. When summed for the 
total stock across the USA GoM, we estimated a stock of 36.4 Tg 
Co, in seagrass sediments from the SE Florida tip to the Rio Grande 


US-Mexico border. 


4.7. Regional GoM Cy,¢ losses 


If the loss estimate is calculated as the “simple” product between 
mean spatial extent and mean global loss (i.e. Duarte, 2002; 
Fourqurean et al., 2012a), the estimated loss for the USA GoM would 
be 17.7 Tg Co... The “simple” estimate differs by almost 20% from 
our “adapted” estimate of 21.69 Tg C,,, found in the upper 20 cm of 
sediment across the USA GoM (Appendix C). We estimate that since 
WWIL, using a series of governmental scientists’ estimates of loss that 
the total spatial seagrass loss in the USA Gulf of Mexico is approxi- 
mately 1,189,244 ha (see Appendix B for extents and references). We 
estimate a USA GoM-wide loss estimate of 21.69 Tg Corg during the 
past 50-75 y (Appendix C), although for Mexico's Yucatan Peninsula 
and northeastern GoM Mexican coast, there are no adequate seagrass 
spatial extent loss data. 


5. Discussion 


Our results demonstrate that sedimentary carbon stocks and fluxes 
vary in accord with the environmental setting of the estuary, species, 
impact type, and subsequent conservation status after restoration. Our 
study also highlights the importance of seagrass restoration as an im- 
portant climate change mitigation strategy that has high potential for 
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Table 2 

Study site sedimentary carbon characteristics for top 20 cm of sediment at the 8 GoM sample sites including organic carbon losses from impacts and gains from seagrass restoration 
measured in and adjacent to selected of the US GoM seagrass restoration sites of Thorhaug (1974, 1985, 2001, 2014) and Thorhaug et al. (1978). Seagrass bed health and protection 
status indicate whether the site has experienced continuous impact exposure over time or if it was subsequently put under conservation status. Impact event age indicates the number 
of years since the impacting event occurred that initially stimulated the seagrass die-off event. 


Seagrass health Impact Total Cog Cog burial rate of Total C,.. Cor, gain per year 
and protection event age loss(Mg_ Restoration restored beds gain (Mg _ since restoration (Mg 
Site status State Sedimentary C,,, stock (Mg ha” 4 (years) ha ') age (years) (gm “y ') ha ') ha !y ') 
Always Impacted Natural Restored 
barren _ barren seagrass seagrass 
Turkey Impacted then Florida 35.61 44.13 50.74 116.06 43 37.26 35 16.828 71.93 2.06 
Point conserved 
27th Street Continually Florida 15.68 31.89 17.73 69.31 57 5.99 32 10.050 37.41 1.17 
impacted 
36th Street Continually Florida 14.58 7.46 53.73 41.18 56 53.36 35 5.971 23.72 0.68 
impacted 
Grove Isle Impacted then Florida 15.61 4.65 34.44 30.77 54 O55 30 4.462 16.13 0.54 
conserved 
Predator Impacted then Texas 6.70 5.77 17.87 13.69 90 7.68 16 2.054 7.92 0.50 
conserved 
San Luis Continually Texas 6.50 0.48 6.89 11.42 55 4.26 1 1.656 0.95 U9. 
Pass impacted 
Shamrock Continually Texas 10.94 3.28 19.54 22.87 65 1.82 iz 3.431 9.60 4.80 
impacted 
Aransas Continually Texas 4.11 4.90 4.46 4.95 10 7.92 1 0.743 0.05 0.05 
Bay impacted 
50 
@ Always barren 4 
0 Natural seagrass 
A |mpacted barren 40 
35 
FS 30 
E $25 
£ oO 
£ om 20 
$ _ 
= 10 
eo 
£ 5 
3 0 
young intermediate old 
seagrass restoration age 
Fig. 5. Eight GoM seagrass restoration sites showing mean (+ SE) organic C,,, Mg 
ha! of young (<2-year-old), intermediate (15—1l6-year-old), and _ older 


0.000 0.003 0.005 0.007 0.010 0.013 0.015 0.018 0.020 
Cog density (g/em3) 


Fig. 4. Vertical depth core profiles of Gulf of Mexico seagrass of mean organic car- 
bon (C,,,) density (g/em*) comparing mean of treatments natural seagrass, always bar- 
ren sites, and previously impacted but remaining barren sediment treatments at 8 USA 
GoM measured seagrass restoration sites. Horizontal bars represent standard errors for 
each depth interval. Vertical bars indicate the midpoint of the depth interval of each 
sample (i.e. 0-5 cm, 5-10 cm, or 10-20 cm). Abbreviations are: AB = always barren 
sites, IB = impacted barren-never-restored or revegetated sites; and NS = natural sea- 
grass. AB C,,, density differs significantly from NS and IB across all sediment depths, 
while IB Corg density is significantly lower than NS at only the shallower sediment 
depths according to mixed model analyses (P < 0.05). 


incorporation into carbon accounting schemes at the national- and re- 
gional-levels. Our regional-scale estimate of naturally-occurring sea- 
grass sedimentary carbon stocks, loss from disturbance, and gains 
through restoration provide a first-time estimate of the Co,.-cycling in 
seagrass across the GoM region. This regional-scale assessment relies 
on a systematic methodology of analysis, which uses a barren con- 
trol to assess baseline C,,, per at each estuary and consistently imple- 
mented restorations at each site. We added the natural barren back- 
ground control as an important element for measuring the amount of 


(38—-43-year-old) restored seagrass sedimentary organic carbon. Older beds had signif- 
icantly higher C,,, stocks than young or intermediate restored sites, from the mixed 
model comparisons (P = 0.0025). 


loss of seagrass beds after impact, and for estimating the actual con- 
tribution of natural seagrass beds and seagrass restoration over back- 
ground sedimentary carbon to the blue carbon budget. 

While Lavery et al. (2013) and Duarte et al. (2013) have noted 
that seagrass blue carbon stocks and metabolism can display signifi- 
cant spatial heterogeneity, we demonstrated that both C,,, losses and 
gains through restoration were similarly influenced by local environ- 
mental variation, by impact type, and by whether the site was sub- 
sequently preserved or continually disturbed after seagrass die-off. 
We showed that restored beds in areas that were initially lost from 
sewage effluents had significantly higher carbon stocks after restora- 
tion than seagrass areas that were lost and restored after dredge and/ 
or fill impacts. We hypothesize this was probably a result of the high 
nitrogen and phosphate inputs associated with the sewage spills re- 
maining in the sediment. Likewise, C,,, stocks were lower in restored 
beds that were continuously disturbed or impacted following restora- 
tion than in beds that were impacted once and then conserved. Our 
ability to delineate this information hinged on our unique study de- 
sign that included two seagrass treatments (natural and restored) and 
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Fig. 6. Mean (+ SE) sedimentary organic carbon stocks of natural vs. impacted barren 
(anthropogenically disturbed but never recolonized or restored) seagrass beds for differ- 
ent disturbance event types as treatments in 8 sampled sites the USA Gulf of Mexico. 
Treatments differed significantly (P < 0.05, depicted with asterisks) within each impact 
event type for dredge and fill, turbidity, vessel grounding, and sewage impact events ac- 
cording to the mixed effect model analyses. 


two controls (historically-barren and presently-barren-previously-im- 
pacted) plus our intensive sampling regime at 8 sites around the USA 
GoM. These results also highlighted differences in background (bar- 
ren) Cor, stocks among estuaries, which suggests that natural baseline 
environmental carbon variation is an important factor in the C-seques- 
tration potential of seagrass beds. For example, the Florida estuaries 
experience largely oceanic C inputs because of the lack of massive 
rivers that drain into those sites, while the Texas sites receive larger 
terrigenous riverine inputs from large watersheds. Likewise, the mag- 
nitude of the combined blue carbon sink differed among sites, with 
some estuaries having extensive adjacent marshes, some with man- 
groves and seagrass, and some with seagrass only, all of which dif- 
ferentially influence C,,, stocks and cycling (Hemminga and Mateo, 
1996; Gonneea et al., 2004; Kennedy et al., 2010). Additional studies 
using stable isotopes are warranted at these sites to explore differences 
in C inputs and cycling among different environmental conditions sur- 
rounding the seagrass beds. 
5.1. Regional differences in C stocks among sites 

The largest C,,. stock in the GoM lies in South Florida in the lower 
subtropical to tropical area between the reef tracks, Florida Keys back 
reef, and the GoM FL shelf and Florida Bay, including the South 
Florida estuaries. In Texas, the highest Co,. stocks are in the Laguna 
Madre, between the Padre Island National Seashore (eastern shore) 
and the massive King Ranch (western shore). For Mexico, precise spa- 
tial extent measurements are being formulated, but are not yet avail- 
able. However, the carbon stock for this portion of the GoM is ex- 
pected to be far smaller than 0.6 Tg. Neither Cuba nor the Mexi- 
can GoM has seagrass restoration attempts to-date. Of the over | mil- 
lion ha of seagrass lost in the GoM, Texas has lost a far greater propor- 
tion of its seagrass than Florida, especially in the northern, industrial- 
ized estuaries. It is clear from Appendix B that Louisiana, Alabama, 
and Mississippi are each close to the “tipping points” of C,,,. loss since 
over 80% or more of their seagrass cover has been lost over the last 
few decades. The State of Florida Department of Environmental Pro- 
tection states losses of over one million ha. We note that if seagrass 
restoration were undertaken for even half of the million ha of sea- 
grass that was previously lost, this could contribute up to an additional 
12.7 Tg Corg to the existing 34.6 Tg of standing stock. 


Such results demonstrate that, in the GoM, seagrass is a major blue 
carbon sink, and that restoration of seagrasses in areas that previously 
hosted seagrass cover in this region can offer a viable climate change 
mitigation strategy. Such activities have great potential for integration 
into carbon accounting schemes to enhance blue carbon conservation 
such as MRV (Monitoring, Reporting, and Verification) and for min- 
imizing C emissions from seagrass loss through mechanisms such as 
REDD + carbon accounting schemes. The utility of seagrass restora- 
tion also ameliorates other ecosystem services, including bolstering 
fisheries, fish nursery function, and stabilizing sediments (Moberg and 
R6nnbick, 2003; Barbier et al., 2011; van Katwijk et al., 2016). This 
study demonstrates that subtropical/tropical restored seagrass can se- 
quester large amounts of carbon within a short time following restora- 
tion (as seen at San Luis Pass and Aransas Bay, Texas), which dif- 
fers from mangroves that take some years before sequestration occurs 
in large amounts (Del Vecchia et al., 2014). These results suggest that 
carbon accumulation in restored beds occurs quickly after restoration, 
likely through a combination of the following processes: direct sedi- 
ment trapping, pumping of photosynthates into the sediments by sea- 
grass plants, and diagenesis. This rapid C,,, sequestration following 
restoration activities highlights the potential of seagrass restoration for 
ameliorating the effects of blue carbon losses in the GoM and else- 
where. 

Florida's large seagrass C,,, sink is probably due to a longer coast 
line, and because only a few short, clastic rivers drain into this part of 
the GoM which can cause severe turbidity and salinity alterations. The 
large seagrass preservation areas set aside since the 1960's (i.e., the 
Everglades National Park, 3 national marine estuaries, Federal EPA 
estuaries, and State preserves) and governmentally-enforced seagrass 
conservation over the last 45 years likely also contributed to the high 
Corg stocks at our Florida samples sites. Florida's program of public 
awareness has created a high value for seagrass among citizens and 
government resource managers, who demanded substantial amounts 
of restoration for seagrass “taking” compared to the other GoM states. 
As a result, Florida has less proportional seagrass loss and more 
ha of restored seagrass than other US GoM states (Thorhaug, 2014; 
Thorhaug et al., in review). Louisiana's far more extensive shoreline 
than Mississippi or Alabama, is dominated by marshes, is well-docu- 
mented as retreating from disturbances, experiences subsidence, and is 
highly vulnerable to sea-level rise and intense storm-events (Handley 
et al., 2007). Louisiana's seagrass is presently limited to an eastern off- 
shore extent behind the eroding Chandeleurs barrier islands in the Bre- 
ton National Wildlife Refuge (Handley et al., 2007). 


5.2. Loss profile for sediment depth of organic carbon 


We hypothesize that a major pathway of the seagrass sediment car- 
bon loss process is characterized by a slow physical diffusion down 
a gradient. The loss of carbon in the upper 10 cm is important to 
the discussion of core length for the “best” carbon loss measure- 
ments by Fourqurean et al. (2012a) suggest that up to 1 m of sed- 
iment may be lost when a seagrass meadow is impacted, but this 
one meter loss estimate remains highly debated. What our current 
study suggests is that rate of egress of carbon in the upper 10 cm 
layer occurred more rapidly than in deeper sediments 10-20 cm af- 
ter the seagrass meadow was killed, and that some the sediment did 
not leave the site. The differences in C,,. among the impacted bar- 
ren core depth profile, the always barren carbon depth profile, and 
the natural seagrass profile indicate that some C,,. loss has occurred 
in the impacted sites at depths upper portions of the sediment profile 
(0-10 cm). However, similarities at the deeper, 10-20 cm core depths 
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among the impacted barren, restored, and natural seagrasses indicate 
that areas that have been previously-disturbed are likely retaining 
some residual carbon from the prior natural seagrass bed. These results 
clearly suggest that residual natural seagrass carbon in 20 cm-deep 
sediments may remain for decades at a site following disturbance. 
This contradicts the Fourqurean et al. (2012a) hypothesis of total loss 
of 1 m of carbon. Variability among sites could explain this discrep- 
ancy, as could differences in carbon loss following various impact 
events. However, there are natural and anthropogenic carbon loss fac- 
tors influencing seagrass C,,, stock and flux across the GoM, and loss 
pathways remain a “black box” of variables without further experi- 
mentation. High energy extreme storm events occur in areas such as 
Florida Bay. Rivers are major sources of CO, efflux (Raymond et al., 
2013), and high-intensity storm events throughout the GoM also de- 
liver large amounts of CO, to coastal marine systems (Bianchi et al., 
2013; Hotchkiss et al., 2015; Raymond and Spencer, 2015). Yet, the 
effect of such sedimentation from storm events on seagrass carbon 
stocks, sequestration, and flux has yet to be estimated. Nor have the 
impacts of sediment reworking and bioturbation on C-cycling caused 
by animals such as Callianassa shrimp (Roche et al., 2016), clams, 
burrowing crabs, rays, and serpulid worms. Sediment microbial activ- 
ity also influences carbon-cycling within natural (Ziegler and Benner, 
1999; Holmer et al., 2004) and restored seagrass beds. Limited data 
from 7. testudinum in Laguna Madre, Texas via phospholipid fatty 
acid concentrations by Jones et al. (2003) suggests that bacterial abun- 
dance is greater in vegetated communities compared to bare areas. The 
microbial community may be responsible for a significant component 
of dissolved organic carbon remineralization in seagrass beds. How- 
ever, future studies that examine the role of bacteria and plankton in 
relation to seagrass Cy, losses from impact and subsequent gains from 
restoration could facilitate a better understanding of seagrass biogeo- 
chemistry in the GoM and elsewhere (e.g., Liu et al. 2017). 

Future research that investigates the rate and magnitude of C,, 
losses across the depth profile of seagrass sediments is needed. In the 
context of seagrass restoration, the more rapidly an impacted site is 
restored, the greater the potential for mitigating C,,, losses, including 
stopping carbon leakage to estuarine waters. Stopping leakage is per- 
tinent given the current debates about how deep Co,, losses penetrate 
into the sediment profile (i.e. Pendleton et al., 2012), and how differ- 
ent environmental settings influence C,,, flux (i.e. Marba et al., 2015). 
It is theoretically possible for the entire meter to be disturbed in very 
high energy sites as in Perth, Australia by Marba et al. (2015). How- 
ever, our study sites occurred in the sheltered, lower energy conditions 
typical of GoM estuarine seagrass meadows, where not all of the top 
20 cm of sediment was lost over the multi-decadal period after the im- 
pacting events. Our restorations dated from events ranging from 1945 
to 2012 during an intense industrialization and urbanization period in 
the GoM, in which impact events were still very apparent in the raw 
sediment profiles. 

The enclosed nature of the GoM's basin with sediment accretion 
throughout soft-bottom estuaries invites comparisons with other areas 
such as the North American Atlantic and Pacific coasts which con- 
tain seagrass at more high energy sites. For example, Lavery et al. 
(2013) also identified variability among Cy, stocks when comparing 
both sites and species at 17 sites within a high energy section of the 
open Australian coastline to the Pacific Ocean, that ranged from tem- 
perate to subtropical. Such results herald the need to incorporate spa- 
tial variation in Co, stocks in larger-scale estimates of seagrass carbon 
stock and sequestration. Spatial interpolators and multivariate model- 
ing could enhance the accuracy of such estimates by providing tools 


that incorporate the range of variation in C 
estimates. 


org Stocks into large-scale 


5.3. Sedimentary carbon loss comparison between species 


Our results suggest that different species can sequester different 
amounts of sedimentary Cy,,. For example, the sediments beneath 7: 
testudinum were more carbon rich than those below H. wrightii which 
corroborates other prior research by Miyajima et al. (2015), Serrano 
et al. (2013) and Lavery et al. (2013) who all demonstrated that C 
stocks varied among Asian and Australian subtropical-Pacific sea- 
grass species. In the GoM, T. testudinum has extensive root and rhi- 
zome systems known to influence C concentrations in seagrass sedi- 
ments (Wetzel and Penhale, 1979; Moriarty et al., 1986). Thus, exper- 
iments that compare carbon translocation by roots to sediment for the 
range of dominant GoM seagrass species are needed for understand- 
ing species-level differences in C,,, sequestration. 


org 


5.4. Impact effects on sequestration 


Low salinity manipulation (diverting water for agriculture, urban- 
ization, etc.) had negative influences on Cy, sequestration at the 
Aransas Bay site. Sewage effluents in sediments were a stimulant to 
seagrass C,,. sequestration, which accompanied high seagrass growth 
after restoration measured at this site (Thorhaug, 1987). Four of our 
sample sites suffered only a single impact period at the time of the 
“killing” of seagrass and are now preserves. Four other sites ex- 
perienced continually high disturbances; these disturbed sites con- 
tained less sedimentary C,,, carbon post-restoration than the four sin- 
gle-event-impact-conserved sites. Obviously, the areas of greatest im- 
pact effects are also where greatest losses of seagrass and its carbon 
have occurred. 


5.5. Rate of seagrass organic carbon sequestration into sediment 


Previous temperate studies of restored vs. natural sequestered sea- 
grass carbon (Greiner et al., 2013; Marba et al., 2015) did not show 
greater sequestration in restored seagrass when compared to natural 
seagrass, while our results indicated that at some sites, restored sea- 
grass C,,. was significantly higher than that of adjacent natural sea- 
grass beds. Our hypothesis is that more vigorous growth of younger 
plants in restored seagrass beds sequesters more sedimentary carbon 
than older natural seagrass meadows. The age of restoration had a sig- 
nificant effect on carbon stocks and C,,, sequestration rates appear 
larger for southern than northern GoM areas, partially due to vigorous 
12-month seagrass growth at temperatures from 21 to 29 °C across the 
sample sites (Thorhaug and Roessler, 1977; Dunton, 1994; Livingston 
et al., 1998). 


6. Conclusions 


Seagrasses provide a multitude of ecosystem services across the 
the United States and Mexico. Our work identifies the GoM as a North 
American seagrass blue carbon hotspot with substantial natural sea- 
grass C,,, sequestration. Due to the already high C,,, stocks in nat- 
ural seagrass beds throughout the region, it also represents a prime 
location for seagrass restoration as a climate change mitigation strat- 
egy. Although, seagrass is abundant across much of the coastal USA, 
where similar calculations of this type could contribute greatly to our 
understanding of the potential for blue carbon accounting and restora- 
tion potential, seagrasses in temperate regions of the US may not 
sequester as much carbon as the subtropical/tropical GoM due to a 
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more limited growing season and a longer history of anthropogenic 
stressors. The IPCC sub-panel on wetlands indicated that seagrass Co,, 
data is needed for the integration of seagrass restoration in REDD + 
carbon mitigation. Carefully chosen sites for well-executed seagrass 
restoration in areas that are likely to be conserved over the long-term 
appear to us as solid mitigation tools for IPCC blue carbon emission 
factor offsets (Thorhaug et al., 2009; van Katwijk et al., 2016). We 
conclude that seagrass must be simultaneously conserved throughout 
the world and be restored in the substantial areas where seagrass has 
previously been lost from anthropogenic activities to avoid further 
sedimentary carbon losses. Both government and private sector sea- 
grass restoration projects will be necessary for mitigating prior sea- 
grass C,,. losses while bolstering other ecosystem services which nat- 
urally accrue from restoration (e.g. fisheries and food-web enhance- 
ment, shoreline resilience, endangered species habitat — for Whooping 
Cranes, Kemp Ridley's turtles and manatee habitat, water clarity, and 
nutrient cycling). 
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Appendix A 


Mean (+ SE) dry bulk density of sediment at different depth below 
the surface in always barren, impacted barren, and natural and restored 
seagrass sediments in the Gulf of Mexico eight restoration sites.! 
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‘tt is estimated that over the last 50 years, seagrass habitat losses ranged 
anywhere from 20% to 100% for most estuaries in the region (Handley et al., 
2007)". In 1992, the Gulf of Mexico Program's Habitat Degradation Committee 
produced a report on the status and trends of emergent and submerged vegetated 
habitats in the Gulf of Mexico (Duke and Kruczynski, 1992). At that time, the total 
seagrass coverage in the shallow estuarine and nearshore waters of the northern 
GoM was estimated to be 1,019,844 ha, but offshore Halophila extent especially 
in federal waters was undefined. Thus, estimates of loss of deep water seagrass are 
very poorly established. 


Appendix B 


Seagrass coverages by states and sites and gains. Data for spatial 
extent from a series of investigators, Eleuterius and Miller (1976); 
Enos and Perkins (1979); Brown and Fisher (1980); Darnell et al. 
(1983); Iverson and Bittaker (1986); Sargent et al. (1995); Pulich 
(1998); Handley et al. (2007); Moncreiff (2007); Fodrie et al. (2010); 
Carter et al. (2011); Lewis (2011); Yarbro and Carlson (2011). Se- 
questration data is from seagrass restorations by Thorhaug (1974); 
Thorhaug (1980a); Thorhaug (1985); Thorhaug (1987); Thorhaug 
(2001); Thorhaug (2014) (8 sites herein), using Fourqurean et al. 
(2012b) data (upper 20 cm) for Florida Bay, southwest Florida, and 
Florida Keys. 


Spatial 
extent Sediment nat- _ Spatially “ad- Total 
(gov't. ural seagrass —justed” totalcar-  C,, ,/area 
mea- measured bon standing Tg Cor, in 
sured) herein (C,,, stock 20 cm in 20 cm from 


ate & sites ectares gha_ area Mg this study 
State & si h Mg ha ! Mg Cy, his stud 


TEXAS 94,389 1,546,949 1.5 
Galveston & 210 4.6 527 
Christmas 
Bays* 
Matagorda. 5840 11.82 69,029 
Bays* 
Sn. Anton. 3238 9.60 31,065 
Bay/Aransas/ 
St Charles* 
Corpus & Red- 9950 11.82 117,668 
fish, Nueces" 
Laguna Madre* 75,151 17.68 1,328,670 
LOUISIANA 1012 37.26 37,707 0.04 
Chandeleurs? 
MISSISSIPPI 750 212 17,900 0.18 
Mississippi 
Sound® 
ALABAMA 202 13.57 2741 0.003 
Coasts/Mo- 
bile’ 
FLORIDA 850,974 34,775,381 34.7 
Panhandle** 15,881 37.4 593,949 
Big Bend*® 248,850 55.66 13,850,991 
Central® 23,790 37.4 889,746 
Southwest® 30,169 32.78 986,526 
Florida Bay*® 58,804 32.78 1,922,890 
Florida Keys/ 346,555 30:78 11,332,349 
GoM shelf® 
Atlantic Shelf* 57,024 37.4 2,132,697 
Turkey Point 37,620 57.05 2,146,221 
Grove Isle 21,533 34.44 741,597 
24-27th St 5374 17.73 95,281 
Julia Tuttle 5374 15.41 83,135 
Total 947,327 36,380,688 Mg Cy. 36.4 Tg 


* Pulich (42). 

» Poirrier and Handley (45). 

® Sediment carbon data from 0.20 from Fourqurean, Kendrick, Collins, Chambers and 
Vanderklift (9). 

© Montcreiff (33). 

® Yarbro and Carlson (34). 

' Carlson and Madley (2006). 


Appendix C 


“Adapted” organic carbon loss estimates calculated from our sea- 
grass organic sedimentary measurement of loss most comparable to 
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each sub-region times spatial loss estimates by government's scien- 
tists. A list of specific zones around the Gulf of Mexico are herein 
used to calculate specific amounts of measured loss by Handley et 
al. (2007), Pulich (1998), Brown and Fisher (1980), Darnell et al. 
(1983), Eleuterius and Miller (1976), Moncreiff (2007), Iverson and 
Bittaker (1986), McNulty et al. (1962), Enos and Perkins (1979) et al.; 
Sargent et al. (1995), and Chardon (1976). Seagrass sedimentary car- 
bon content data of natural seagrass beds, minus the sum of (previ- 
ously impacted seagrass bed minus always barren sediment) is from 
Thorhaug (1980b), Thorhaug (1985), Thorhaug (1987), Thorhaug 
(2001), Thorhaug (2014), Thorhaug et al. (in review) (8 sites herein), 
using Fourqurean et al. (2012b) carbon data (upper 20 cm) for Florida 
Bay and surrounds. 


Loss 
mea- 
sured 
herein Spatially Total 
Loss spa-___rele- adjusted Corg 
tial vant total loss loss/Tg 
(govt. Corg in stand- Corg in 
mea- Mg ing stock 20 cm 
sured) ha ‘in (Mg Corg depth 
in in in 20 cm/ per re- 
State & site hectares 20 cm region) gion Source 
TEXAS 8,203,618 8.26 
Galveston 110,365 4.26 470,154 Pulich (1998) 
Christmas 
Bays 
Matagorda 184,285 11.82 2,178,249 Pulich (1998) 
Bays 
Sn Antonio/ — 91,983 9.60 883,037 Pulich (1998) 
Aransas/St 
Charles 
Corpus/Red- 78,931 11.82 932,964 Pulich (1998) 
fish 
Laguna 211,494 17.68 3,739,214 Pulich (1998) 
Madre 
LOUISIANA = 5366 13.57 72,817 0.07 Poirrier and 
Tles de Chan- Handley (2007) 
deleurs 
MISSISSIPPI 6127 21.2 129,892 0.13 Duke and 
Mississippi Kruezynski 
Sound (1992); Handley 
et al. (2007), 
Moncreiff 
(2007), 
Eleuterius and 
Miller (1976), 
Carter et al. 
(2011) 
ALABAMA — 930 4.26 4054 0.004 Iverson and 
Coasts & Bittaker (1986); 
Mobile Fodrie et al. 
(2010); Heck et 
al. (2015) 
FLORIDA 13,389,440 13.39 
North Pan- 5335 13.57 72,395 Handley et al. 
handle (2007), Lewis 
(2011) 
Big Bend 327,939 22.47 7,368,767 Mattson (1999), 
Duke and 
Kruezynski 
(1992) 
Iverson & Bit- 
taker 
Central 22,903 37.26 851,991 Handley et al. 
(2007), Duke and 
Kruezynski 
(1992) 
Iverson & Bit- 
taker 
Southwest 55,881 32.7° 1,827,309 Handley et al. 
(2007) 


Florida Bay 27,000 32.7" 882,900 Robblee et al. 
(1991), Handley 
et al. (2007), 
Fourqurean et al. 
(2012b) 
Florida Keys 29,749 37.26 1,108,448 Handley et al. 
(2007), 
Atlantic shelf 5703 37.26 212,440 Handley et al. 
(2007) 
Card Sound 10,000 57.12 571,200 Handley et al. 
& south Bis- (2007), 
cayne Bay (Thorhaug and 
Roessler, 1977); 
Thorhaug 
(1980a) 
Central Bis- 5755 22.47 129,315 Harlem (1979), 
cayne Bay Handley et al. 
(2007), 
Thorhaug 
(1980b); 
Thorhaug 
(1987), Chardon 
(1976), 
Lower half 4749 50.44 239,539 Thorhaug et al. 
North Bis- (1979), Chardon 
cayne Bay. (1976), 
Upper half N. 4749 3.14 14,912 Chardon (1976), 
B.B Thorhaug (1987) 
Total loss 1,189,244 21,689,597 21.69 Tg 


** Data from Fourqurean et al. (2012a) is used herein in Florida Bay and Southwest 
Florida for the upper 20 cm of the total 1 m of the core (163.5 Mg for 1 m) to make data 
consistent with the loss sampling unit of 20 cm. 
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